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We present the results of a search for the production of an excited state of the electron, e* , in 
proton-antiproton collisions at yfs = 1.96 TeV. The data were collected with the DO experiment at 
the Fermilab Tevatron Collider and correspond to an integrated luminosity of approximately 1 fb -1 . 
We search for e* in the process pp — > e*e, with the e* subsequently decaying to an electron plus 
photon. No excess above the standard model background is observed. Interpreting our data in the 
context of a model that describes e* production by four-fermion contact interactions and e* decay 
via electroweak processes, we set 95% C.L. upper limits on the production cross section ranging 
from 8.9 fb to 27 fb, depending on the mass of the excited electron. Choosing the scale for contact 
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interactions to be A = 1 TeV, excited electron masses below 756 GeV are excluded at the 95% C.L. 
PACS numbers: 12.60.Rc, 14.60.Hi, 12.60.-i, 13.85.Rm 



An open question in particle physics is the cause of the 
observed mass hierarchy of the quark and lepton SU(2) 
doublets in the standard model (SM). One proposed ex- 
planation for the three generations is a compositencss 
model |l| of the known leptons and quarks. According 
to this approach, a quark or a lepton is a bound state of 
three fermions or of a fermion and a boson Due to 
the underlying substructure, compositeness models im- 
ply a large spectrum of excited states. The coupling of 
excited fermions to ordinary quarks and leptons, result- 
ing from novel strong interactions, can be described by 
contact interactions (CI) with the effective four-fermion 
Lagrangian [3J] 
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where A is the compositeness scale and is the fermion 
current 

ip. = vl hi^fh + v'l flinfl + vl flltifh 

+ h.c. + (£-> R). 

The SM and excited fermions are denoted by / and /*, 
respectively; g 2 is chosen to be 47r, the 77 factors for the 
left-handed currents are conventionally set to one, and 
the right-handed currents are set to zero. 

Gauge mediated transitions between ordinary and ex- 
cited fermions can be described by the effective La- 
grangian [H 
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where G^ v , W^ Vl and are the field strength tensors of 
the gluon, the SU(2) and U(l) gauge fields, respectively, 
and f si f and /' are parameters of order one. 

For the present analysis, we consider single produc- 
tion of an excited electron e* in association with an elec- 
tron via four-fermion contact interactions, with the sub- 
sequent electroweak decay of the e* into an electron and a 
photon (Fig. ma)). This decay mode leads to the fully re- 
constructable and almost background-free final state eej. 
With the data considered herein, collected with the DO 
detector at the Fermilab Tevatron Collider in pp collisions 
at y/s = 1.96 TeV, the largest expected SM background is 
from the Drell-Yan (DY) process pp — > Z/-f* — > e + e~('y), 
with the final state photon radiated by either a parton in 
the initial state or from one of the final state electrons. 
This background can be strongly suppressed by the appli- 
cation of suitable selection criteria. Other backgrounds 
are small. 

Previous searches have found no evidence for the pro- 
duction of excited electrons, e.g. at the CERN LEP e + e~ 



H and the DESY HERA ep H colliders, in the con- 
text of models where the production of excited electrons 
proceeds via gauge interactions; however, the reach has 
been limited by the available center-of-mass energy to 
m e* ^ 300 GeV. Searches for quark-lepton composite- 
ness via deviations from the Drell-Yan cross section at 
the Tevatron have excluded values of A of up to « 6 TeV 
depending on the chirality @. The present analysis is 
complementary to those results in the sense that an ex- 
clusive channel and different couplings (j) factors) are 
probed. The CDF collaboration has recently presented 
results [7| for the production of excited electrons which 
will be discussed later. 

(a) (b) 
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FIG. 1: (a) Four-fermion contact interaction qq — > e*e, and 
electroweak decay e* — > ery. (b) Relative contributions of de- 
cays via contact interactions and via electroweak interactions 
(GM) as a function of m e « /A. 

For the simulation of the signal the pythia event gen- 
erator [§[ is used, following the model of Ref. ||. The 
branching fraction for the decay e* — > ej normalized 
to all gauge particle decay modes is 30% for masses 
above 300 GeV; for smaller e* masses it increases up 
to 73% at m e * = 100 GeV. Decays via contact inter- 
actions, not implemented in pythia, contribute between 
a few percent of all decays for A ^> m e » and 92% for 
A = m e . [| (see Fig.QJb)). This is taken into account for 
the signal expectation. The leading order cross section 
calculated with pythia is corrected to next-to-next-to- 
leading order (NNLO) using Ref. 9]; the corresponding 
correction factor varies between 1.37 and 1.42, depend- 
ing on the invariant mass of the electron and the ex- 
cited electron. The total width is greater than 1 GeV for 
100 GeV < m e . < 1000 GeV, thus lifetime effects can be 
neglected. For the values of m e * and A studied here, the 
total width is always less than 10% of m e * Q. 

The dominant SM background process at all stages 
of the selection is DY production of e + e~ pairs. This 
background, as well as diboson (WW, WZ, ZZ) produc- 
tion, is simulated with the pythia Monte Carlo (MC) 
program. The DY expectation (as well as W — ► ev) 
is corrected using the NNLO calculation from Ref. 0. 
For diboson production, the next-to-leading order cross 



sections from Ref. [Tcj are used. Contributions from tt 
[ll| and W boson production are found to be negligible. 
Monte Carlo events, both for SM and signal, are passed 
through a detector simulation based on the GEANT [jjj 
package and reconstructed using the same reconstruction 
program as the data. The CTEQ6L1 parton distribution 
functions (PDFs) [l3| are used for the generation of all 
MC samples. 

The analysis is based on the data collected with 
the DO detector [T3] between August 2002 and Febru- 
ary 2006, corresponding to an integrated luminosity of 
1.01 ± 0.06 fb -1 . The DO detector includes a central 
tracking system, which comprises a silicon microstrip 
tracker and a central fiber tracker, both located within 
a 2 T superconducting solenoidal magnet, and optimized 
for tracking and vertexing capability at pseudorapidities 1 
1 77 1 < 2.5. Three liquid argon and uranium calorime- 
ters provide coverage out to |ry| w 4.2: a central sec- 
tion (CC) covering \rj\ < 1.1, and two end calorime- 
ters (EC). The electromagnetic section of the calorimeter 
has four longitudinal layers and transverse segmentation 
of 0.1 x 0.1 in 77 — cj) space, except in the third layer, 
where it is 0.05 x 0.05. A muon system resides beyond 
the calorimetry, and consists of layers of tracking detec- 
tors and scintillation trigger counters before and after 
1.8 T iron toroids. Luminosity is measured using scin- 
tillator arrays located in front of the EC cryostats, cov- 
ering 2.7 < \rj\ < 4.4. A three-level trigger system uses 
information from tracking, calorimetry, and muon sys- 
tems to reduce the pp bunch crossing rate of 1.5 MHz to 
~ 100 Hz, which is written to tape. 

Efficiencies for electron and photon identification and 
track reconstruction are determined from the simulation. 
To verify the simulation and to estimate systematic un- 
certainties, the efficiencies are also calculated from data 
samples, using Z — > e + e~ candidate events and other 
dilepton events for electrons and tracks. Small differences 
between the efficiency determinations from data and sim- 
ulation are corrected in the simulation. We assume that 
the different response for electrons and photons in the 
calorimeter is properly modeled by the simulation. The 
transverse (with respect to the beam axis) momentum 
resolution of the central tracker and the energy resolu- 
tion of the electromagnetic calorimeter are tuned in the 
simulation to reproduce the resolutions observed in the 
data using Z — > 11 (£ = e, p) events. 

The process pp — > e*e with e* — * leads to a final 
state with two highly energetic isolated electrons and a 
photon. First, the two electrons are identified as clus- 
ters of calorimeter energy with characteristic longitudi- 
nal and transverse shower shapes and at least 90% of 
the energy deposited in the electromagnetic section of 
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1 The pseudorapidity r\ is defined as r\ = — In [tan(t9/2)]. We use 
the polar angle 9 relative to the proton beam direction, and <j> is 
the azimuthal angle, all measured with respect to the geometric 
center of the detector. 
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FIG. 2: Invariant dielectron mass distribution in the dielec- 
tron data sample compared to the SM expectation (a) for 
events with both electrons reconstructed in the CC and (b) 
for events in the CC/EC topology, for data (points with sta- 
tistical uncertainties) and SM backgrounds (DY, diboson, tt, 
and multijet production). 



the calorimeter. Two electrons, with transverse energies 
E T > 25 GeV and E T > 15 GeV, are required. Both 
electrons are matched to tracks in the central tracking 
system, and we distinguish between CC (\r]\ < 1.1 with 
respect to the detector center) and EC (1.5 < |?7| < 2.5) 
electrons. Events with the two electrons in opposite EC 
are rejected in order to suppress multijet background. 
The signal is expected to produce isolated electrons, 
therefore both electrons need to fulfill I = (£ ; tot(0.4) — 
£em(0.2))/£U(0.2) < 0.2, where E tot (0A) and E em (0.2) 
denote the energies deposited in the calorimeter and de- 
posited in only its electromagnetic section, respectively, 
in cones of size ATI = \/(Arj) 2 + (A0) 2 = 0.4 and 0.2. 
The electrons are required to be separated by A1Z > 0.4. 

The events were collected with trigger conditions re- 
quiring one or two electrons detected in the calorimeters, 
with varying conditions depending on the Et thresholds, 
the shower shape, tracks in the central tracking system, 
and number of electrons. The overall trigger efficiency is 
determined from independent data samples and is consis- 
tent with 100% for the signal after application of all se- 
lection criteria. The selected dielectron sample contains 
62930 events, whereas 61900 ± 5700 events are expected 
from SM processes. The invariant dielectron mass dis- 
tributions for CC/CC and CC/EC topologies are shown 
in Fig. [2j The largest SM contribution is DY produc- 
tion of e + e~ pairs, followed by multijet production with 
misidentified electrons. The multijet background is es- 
timated using a data sample where at least one of the 
electron candidates fails the shower shape requirements. 
This sample is then corrected as a function of Et and 
?y of the misidentified electrons in order to account for 
different misidentification rates in the CC and EC and 
the different trigger efficiency for misidentified electrons. 

Next, a photon is identified in the event as an iso- 
lated cluster of calorimeter energy with at least 97% of 
its energy deposited in the electromagnetic section of 
the calorimeter (CC or EC). The isolation condition is 
I < 0.07. The photon candidate Et must be larger than 
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FIG. 3: For the eey sample, (a) the photon Et distribu- 
tion, (b) the distribution of the e 2 7 invariant mass com- 
pared with the SM expectation and a possible e* signal 
for m e * = 100 GeV; (c) the e*y invariant mass for the ey 
combination closest to m e » = 800 GeV. In (d) the separa- 
tion A_R(e 2 ,7) is shown after the cut on the invariant mass 
m(e2,7) > 90 GeV for m e * = 100 GeV. The signal corre- 
sponds to A = 2 TeV, 1 TeV, 1 TeV and 4 TeV in (a), (b), (c) 
and (d), respectively. All uncertainties are statistical only. 



15 GeV, no track is allowed to be matched to the pho- 
ton candidate in r\ and <\> with a \ 2 probability of greater 
than 0.1%, and the sum of the transverse momenta of 
tracks within a hollow cone defined by 0.05 < Ai? < 0.4 
around the photon direction has to be below 2 GeV to 
further ensure isolation. Additional shower shape criteria 
are imposed to increase the photon purity. The photon 
candidate is required to be separated from the electron 
candidates in the event by A1Z > 0.4. 

After this selection, 239 ± 36 events are expected from 
SM processes. Of these, 226 events are due to DY 
— > e + e~ with a genuine high Et photon, followed by 
7 events from DY — > e + e~4- jets, where a jet is misiden- 
tified as a photon. The absolute rate of the latter pro- 
cess has been determined from a data sample enriched 
in "fake" photons, applying the rate for such objects to 
be misidentified as photons as a function of Et, and sub- 
tracting the true photon contribution [lj| ■ Finally, about 
4 and 2 events are expected from multijet and diboson 
production, respectively. In the data, 259 events are se- 
lected, compatible with the SM prediction. The photon 
Et distributions for the data and SM background are 
shown in Fig. El^a). 

Additional selection criteria depending on the hypo- 
thetical e* mass are applied to reduce the remaining SM 



background. The following criteria have all been opti- 
mized to achieve the best expected upper limit on the 
production cross section. The e* invariant mass can be 
reconstructed from one of the electrons and the photon. 
For m e * < 300 GeV, the lower Et electron (e2) is pre- 
dominantly the decay electron. Alternatively, of the two 
possibilities to reconstruct the ej invariant mass, the 
value closest to the e* candidate mass is chosen. The 
latter has been done for m e « > 300 GeV. Example mass 
distributions for the two chosen options to reconstruct 
the e* candidate mass are shown in Figs. (3]Jb) and (c). 
The alternatives of single-sided mass cuts and a mass 
window are considered, leading to single-sided cuts for 
all values of m e * . Rejecting events with both electrons 
or the photon in the EC leads to a slightly better sensitiv- 
ity; since for high values of m e * the SM backgrounds are 
extremely small, we have not applied these selection cri- 
teria for m e * > 400 GeV, in order to keep the search gen- 
eral beyond the specific model considered here. Finally, 
the separation A7\L(e2,7) between the lower Et electron 
and the photon allows discrimination between signal and 
background for m e » < 200 GeV. This is illustrated in 
Fig. E^d) for m e * =100 GeV. All mass-dependent selec- 
tion criteria are summarized in Table HI 



TABLE I: Mass-dependent selection criteria. The second and 
the third columns show the lower mass cuts. The next two 
columns show if events with both electrons or the photon in 
the EC are kept, respectively, and in the last column the 
upper value for the separation between the second electron 
and the photon is given. 





m(e 2 ,7) 


m ( e7 )ClOBCSt 


EC/EC EC 


Aft(e 2 ,7) 


[GeV] 


[GeV] 


[GeV] 


e 


7 


100 


> 90 




no 


no 


< 1.8 


200 


> 165 




no 


no 


< 3.3 


300 




> 285 


no 


no 


any 


400 




> 370 


yes 


yes 


any 


500 




> 445 


yes 


yes 


any 


600 




> 515 


yes 


yes 


any 


700 




> 600 


yes 


yes 


any 


800 




> 705 


yes 


yes 


any 


900 




> 800 


yes 


yes 


any 


1000 




> 900 


yes 


yes 


any 



The final selection efficiency varies from 13% (m e » = 
100 GeV) up to « 33% for higher values of m e * . In 
the data we find one event each for the m e * = 200 GeV 
mass hypothesis and for the m e * = 300 GeV mass hy- 
pothesis, respectively, and no events for other values of 
m e * , compatible with the SM expectation. This result is 
summarized in Table HTl 

The dominant systematic uncertainties are as follows. 
The uncertainty on the SM cross sections is dominated 
by the DY process and the uncertainty from the choice 
of PDF and renormalization and factorization scales [(3 
- 10)%]. Electron reconstruction and identification have 
an uncertainty of 2.5% per electron, and a (1 - 4)% uncer- 
tainty is assigned to the photon identification, depending 
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TABLE II: For different values of the e* mass hypothesis, the 
number of selected data events, the SM expectation including 
statistical and systematic uncertainties, and the signal efH- 



5" 10 5 



ciency. 


m e . 
[GeV] 


Data 


OJ-VJ. CAUCLldllUll 




Signal eff. 

[%] 


100 





0.33 ±0.09 ±0.03 




13.2±0.6±1.3 


200 


1 


0.52±0.16±0.05 




16.5 ±0.6 ±1.6 


300 


1 


0.32 ±0.12 ±0.03 




22.2±0.7±2.2 


400 





0.26 ±0.11 ±0.03 




28.3 ±0.8 ±2.8 


500 





0.12 ±0.08 ±0.01 




31.5±1.0±3.1 


600 





(0.57 ±0.54 ±0.06) xl0~ 


i 


32.3 ±0.9 ±3.2 


700 





(0.82 ±0.37 ±0.09) xl0~ 


-3 


34.3 ± 1.1 ±3.4 


800 





(0.48 ±0.28 ±0.06) xl0~ 


-3 


32.2 ±0.8 ±3.2 


900 





(0.17±0.17±0.02) xl0~ 


-3 


33.2 ±0.8 ±3.3 


1000 





(0.17±0.17±0.03) xlO" 


-3 


33.3 ±0.9 ±3.3 



on r\ and Er- The trigger efficiency is 100^3%. The in- 
tegrated luminosity is known to a precision of 6.1% [l6| . 
The uncertainty on the number of background events due 
to jets misidentified as photons is estimated to be 60% 
of this particular background, from differences between 
the expectation from the simulation and the independent 
measurement from the data. A 25% uncertainty is de- 
termined on the multijet background by comparing the 
resulting multijet background estimate when using dif- 
ferent criteria to select the multijet background sample; 
after all selections, the multijet background is negligible. 
The uncertainty on the signal cross section is estimated 
to be 10%, consisting of PDF uncertainties and missing 
higher order corrections. 

Since the observed number of events is in agreement 
with the SM expectation, we set 95% C.L. limits on the 
e* production cross section times the branching fraction 
into e7- A Bayesian technique [i~7| is used, taking into 
account all uncertainties and treating them as symmetric 
for simplicity. The resulting limit as a function of m e « is 
shown in Fig. [4] together with predictions of the contact 
interaction model for different choices of the scale A. For 
A = 1 TeV (A = m e .), masses below 756 GeV (796 GeV) 
are excluded. In Fig. [SJ the excluded region in terms of 
A and m e » is shown. 

The CDF collaboration has recently searched [7J] for 
the production of excited electrons using a data sample 
corresponding to an integrated luminosity of 202 pb _1 , 
but the CDF mass limit of m e » > 879 GeV at the 95% 
C.L. for A = m e * cannot be directly compared to ours for 
two reasons. The e* production cross section calculated 
with the version of PYTHIA used in Ref. [7j is a factor of 
two higher than in subsequent versions corrected by the 
PYTHIA authors. Furthermore, in Ref. @, it is assumed 
that decays via contact interactions can be neglected, 
while in our analysis such decays are taken into account 
in the calculation of the branching fraction e* — > e*y, 
following Ref. [3|. Omitting contact interaction decays, 
we would obtain a limit of m e « > 946 GeV for A = m P . 
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FIG. 4: The measured and expected limits on cross section 
times branching fraction, compared to the contact interaction 
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oretical uncertainty of the model prediction is indicated by 
shaded bands. 
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at the 95% C.L. 2 . 

In summary, we have searched for the production of 
excited electrons in the process pp — > e*e with e* — > 
ej, using about 1 fb _1 of data collected with the DO 
detector. We find zero or one event in the data depending 
on the mass of the hypothetical e*, compatible with the 
SM expectation. We set limits on the production cross 
section times branching fraction as a function of m e * . For 
a scale parameter A = 1 TeV, masses below 756 GeV are 
excluded, representing the most stringent limit to date. 

We thank the staffs at Fermilab and collaborating in- 
stitutions, and acknowledge support from the DOE and 
NSF (USA); CEA and CNRS/IN2P3 (France); FASI, 
Rosatom and RFBR (Russia); CAPES, CNPq, FAPERJ, 
FAPESP and FUNDUNESP (Brazil); DAE and DST 



2 Multiplying the theoretical prediction in addition by a factor of 
two, the mass limit would increase to 989 GeV. 



8 



(India); Colciencias (Colombia); CONACyT (Mexico); 
KRF and KOSEF (Korea); CONICET and UBACyT 
(Argentina); FOM (The Netherlands); Science and Tech- 
nology Facilities Council (United Kingdom); MSMT and 
GACR (Czech Republic); CRC Program, CFI, NSERC 



and WestGrid Project (Canada); BMBF and DFG (Ger- 
many); SFI (Ireland); The Swedish Research Council 
(Sweden); CAS and CNSF (China); and the Alexander 
von Humboldt Foundation. 



[a] Visitor from Augustana College, Sioux Falls, SD, USA. 

[b] Visitor from The University of Liverpool, Liverpool, UK. 

[c] Visitor from ICN-UNAM, Mexico City, Mexico. 

[d] Visitor from II. Physikalisches Institut, Georg- August- 
University, Gottingen, Germany. 

[e] Visitor from Helsinki Institute of Physics, Helsinki, Fin- 
land. 

[f] Fermilab International Fellow. 
[J] Deceased. 

[1] H. Terazawa, M. Yasue, K. Akama and M. Hayashi, 
Phys. Lett. B 112, 387 (1982); F.M. Renard, Nuovo 
Cimento A77, 1 (1983); A. De Rujula, L. Ma- 
iani and R. Petronzio, Phys. Lett. B 140, 253 
(1984); E.J. Eichten, K.D. Lane and M.E. Peskin, 
Phys. Rev. Lett. 50, 811 (1983). 

[2] H. Terazawa, Y. Chikashige and K. Akama, Phys. Rev. D 
15, 480 (1977); Y. Ne'eman, Phys. Lett. B 82, 69 (1979). 

[3] U. Baur, M. Spira and P.M. Zerwas, Phys. Rev. D 42, 
815 (1990). 

[4] C. Adloff et al. (HI Collaboration), Phys. Lett. B 548, 
35 (2002); S. Chekanov et al. (ZEUS Collaboration), 
Phys. Lett. B 549, 32 (2002). 

[5] G. Abbiendi et al. (OPAL Collaboration), Phys. Lett. B 
544, 57 (2002); P. Achard et al. (L3 Collaboration), 
Phys. Lett. B 568, 23 (2003); P. Abreu et al. (DELPHI 
Collaboration), Eur. Phys. J. C 8, 41 (1999); R. Barate 
et al. (ALEPH Collaboration), Eur. Phys. J. C 4, 571 



[6 

[<: 

[§ 
[9 

[10 

[11 

[12 

[13 

[14; 
[is; 

[16 
[17 



(1998). 

B. Abbott et al. (DO Collaboration), Phys. Rev. Lett. 82, 
4769 (1999). 

D. Acosta et al. (CDF Collaboration), 
Phys. Rev. Lett. 94, 101802 (2005). 

T. Sjostrand et al., Comput. Phys. Commun. 135, 238 
(2001). Pythia version 6.3 is used throughout. 
R. Hamberg, W.L. van Neerven and T. Mat- 
suura, Nucl. Phys. B359, 343 (1991); Erratum 
Nucl. Phys. B644, 403 (2002). 

J.M. Campbell and R.K. Ellis, Phys. Rev. D 60, 
113006 (1999); J.M. Campbell and R.K. Ellis, 
http://mcfm.fnal.gov/ 



N. Kidonakis and R. Vogt, Int. J. Mod. Phys. A 20, 3171 
(2005). 

R. Brun and F. Carminati, CERN Program Library Long 
Writeup W5013, 1994 (unpublished). 
J. Pumplin et al, J. High Energy Phys. 0207, 012 (2002); 
D. Stump et al, J. High Energy Phys. 0310, 046 (2003). 
V. Abazov et al. (DO Collaboration), Nucl. In- 
strum. Methods Phys. Res., Sect. A 565, 463 (2006). 
V. Abazov et al. (DO Collaboration), Phys. Lett. B 653, 
378 (2007). 

T. Andeen et al, Fermilab-TM-2365 (2007). 

I. Bertram et al. (DO Collaboration), Fermilab- TM-2104 

(2000). 



